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The  Problem  of  Analog  Fault  Anal  vs  i 

Analog  fault  analysis  is  a  method  of  finding  a  description  (either 

in  terms  of  the  physical  analog  system  or  its  model)  of  the  way  in  which 

an  analog  system  has  failed.  A  failure  (or  fault)  may  manifest  itself  in 

many  different  ways,  depending  on  the  nature  of  the  system.  In  the  present 

study  we  consider  the  system  to  be  a  linear,  analog  electronic  circuit, 

with  at  least  two  accessible  terminals.  We  treat  the  accessible  terminals 

as  a  multi-port  network.  The  behavior  of  the  network  can  be  determined  by 

measuring  each  of  the  several  independent  (complex)  multi-port  parameters 

(transfer  functions),  possibly  at  several  frequencies.  Thus,  if  m/2 

independent,  complex  measurements  are  made,  the  network  performance  is 

given  in  terms  of  m  real  numbers,  each  of  which  is  a  function  of  the 

parameters  of  the  network.  We  denote  the  measurement  set  by  the  real 
*  A  * 

vector  M  -  col  (M. ).  Similarly,  we  denote  the  corresponding  set  of  multi - 
port  parameters  computed  from  the  network  model  by  the  real  vector 
M  a  col  (M. ) . 

We  will  assume  there  are  n  real  network  parameters,  each  of  which  is 
either  a  resistance,  capacitance,  inductance,  or  controlled  source  gain. 

We  denote  the  parameter  set  by  the  real  vector  xj  col  (x.). 

The  fault  analysis  problem  may  be  stated  as  follows:  given  m  real 
measurements  and  the  network  model  containing  n  real  parameters,  find  the 
values  of  the  n  real  parameters  which  produce  the  m  real  measurements . 

This  requires  solving  the  system  of  in  non-linear  equations  in  n  unknowns 

M  ( x )  s  M  ( 1 ) 
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Currently  available  analog  automatic  test  equipment  ( AT C )  allows  the 
preprogramming  of  a  sequence  of  test  signals  and  measurements  [1],  The 
software  interpretation  of  these  measurements  is  mostly  limited  to  a 
functional  test  of  the  component  or  modulo  tested.  To  achieve  greater 
resolution  of  the  fault  with  currently  available  ATL  (i.e.,  to  find  the 
faulty  component  or  integrated  circuit  in  the  module),  requires  access  to 
more  terminals  (test  points)  and  the  making  of  a  more  complex  sequence  of 
measurements . 

Often  in  practice  the  number  of  terminals  available  for  measurement 
is  limited.  For  example,  in  troubleshooting  existing  equipment,  available 
test  points  are  fixed  by  the  design.  In  the  design  of  new  equipment,  it 

is  desirable  to  keep  the  number  of  tost  points  as  small  as  possible  to 

avoid  an  excessive  number  of  pins  in  the  connectors. 

The  present  study  was  directed  at  improving  the  reliability  of  high 
resolution  (down  to  the  component  or  node  if  necessary)  analog  fault 
analysis  from  measurements  at  the  accessible  terminals.  The  measurements 
In  this  study  require  applying  sinusoidal  sources  and  finding  the  magnitude 
and  phase  relationships  of  the  steady  state  responses,  i.e.,  the  measure¬ 
ments  are  a  set  of  complex  multi-port  parameters  at  the  accessible  terminals. 
The  fault  is  isolated  by  solving  (1)  for  the  n-dimensional  parameter  vector 

x.  Recent  work  in  analog  fault  analysis  has  been  directed  toward  finding 

an  efficient  way  to  solve  these  equations. 

1 . 1  Summary  of  Previous  Work 

In  solving  (1),  the  number  of  network  parameters  may  be  greater 
than,  equal  to,  or  less  than  the  number  of  measurements  m.  If  m  <  n 
it  Is  not  possible  to  solve  for  the  parameter  set  x,  since  in  this 
case  (1)  usually  has  infinitely  many  solutions.  Ransom  and  Sacks  [2] 
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approached  this  problem  by  findinp  the  solution  to  (1)  which  minimizes 
the  norm  1 1  x  -  x^  1 1 ,  where  x^  is  the  nominal  parameter  set.  The 
difficulty  with  this  approach  is  that  the  solution  assumes,  rouqhly 
speaking,  that  the  most  likely  state  of  the  clement  values  in  the  net¬ 
work  is  the  one  which  causes  the  smallest  drift  from  the  nominal  values 
consistant  with  the  measurements.  This  assumption  excludes  the  possi¬ 
bility  of  catastrophic  faults  (open  and  short  circuits). 

If  m  £  n,  then  the  non-linear  equations  given  by  (1)  are  often 
solved  using  optimization  techniques  by  finding  min  1 1  M  (x_)  -  H  |  | , 
where  sometimes  the  norm  may  be  reduced  to  zero  if  m  =  n.  This  approach 
was  taken  by  Chen  and  Sacks  [3],  using  an  efficient  algorithm  to 
evaluate  M  (x).  In  order  to  accommodate  catastrophic  faults  (open  and 
short  circuits)  by  this  method,  each  element  in  the  parameter  space 
must  be  searched  over  the  range  (0,  ~) ,  which  is  not  practical  because 
of  the  prohibitively  large  number  of  evaluations  of  the  function  M  (x) . 

In  this  project  we  continue  to  seek  a  way  to  overcome  the  diffi¬ 
culties  of  the  two  methods  described  above,  viz.  that  catastrophic 
faults  arc  overlooked  altogether  or  that  they  are  found  at  the  cost  of 
very  lengthy  computation.  The  approach  described  in  this  report  is  to 
modify  the  parameter  set  x  to  include  the  set  of  all  single  short 
circuits  involving  inaccessible  nodes.  The  remaining  parameters  are 
branch  admittances  or  controlled  source  gains.  Searching  over  short 
circuits  eliminates  the  need  to  search  over  very  large  values  of  the 
branch  admittances.  It  is  also  unnecessary  to  search  over  very  large 
values  of  controlled  source  gains,  since  this  mode  of  failure  is 
extremely  unlikely.  Johnson  [4],  [5]  has  described  a  very  efficient 
algorithm  for  evaluating  the  function  M  (x.)  in  the  presence  of  a  short 
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circuit.  The  efficiency  results  from  the  fact  that  the  network 
change  caused  by  the  short  circuit  results  in  an  alteration  of  the 
(inaccessible)  nodal  -admittance  matrix  by  a  matrix  of  unit  rank. 

As  a  result,  the  response  M  (x)  of  the  faulted  network  can  be 
computed  (in  terms  of  the  nominal  response)  using  many  fewer  multi¬ 
plications  than  arc  required  for  an  original  analysis.  Sparse  matrix 
techniques  were  also  used  to  reduce  computation  time.  The  efficiency 
of  this  algorithm  makes  it  practical  to  exhaustively  search  for  a 
minimum  of  | |  M  (x)  -  M  ||  over  all  possible  single  short  circuits 
Involving  an  inaccessible  node.  The  test  examples  described  by 
Johnson  [4],  [5]  show  that  short  circuits  can  be  reliably  located 
by  the  algorithm,  even  though  measurements  were  made  only  at  one 
frequency,  and  parameter  drift  was  not  taken  into  account. 

Performing  the  measurements  (and  analysis)  at  more  than  one 
frequency  increases  the  number  of  measurements  m  without  changing 
the  number  of  network  parameters  n.  This  may  improve  the  reliability 
of  the  fault  analysis  if  all  the  measurements  are  independent.  The 
problem  of  independent  measurements  is  discussed  in  some  detail  by 
Sen  and  Saeks  [6]  -  [i>],  who  unfortunately  stop  short  of  describing 
exactly  how  to  choose  the  measurements  and  frequencies. 

2 «  New  Results  in  Analog  Fault  Analysis 

The  study  undertaken  by  the  author  was  directed  toward  improving  the 
reliability  of  analog  fault  analysis  (isolating  short  circuits,  open 
circuits  and  parameter  changes)  while  keeping  the  algorithm  as  fast  and 
efficient  as  possible.  Specifically  two  separate  studies  were  done: 
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(a)  Introducing  parameter  changes  Into  the  current  algorithm. 

The  current  algorithm  searches  for  a  minimum  of  | |  M  (x)  -  M  | | 
over  short  circuits  only.  By  searching  over  the  network  parameters 
(resistances,  capacitances,  inductances  and  controlled  source 
gains)  as  well,  the  algorithm  should  be  improved,  because  of  the 
expanded  description  of  the  fault.  In  addition  it  was  hoped  that 
the  algorithm  would  be  made  more  reliable  (i.e.,  it  was  hoped 
that  the  algorithm  would  not  indicate  faults  which  were  not 
present,  but  instead  would  find  only  those  that  were  present). 

An  efficient  algorithm  for  calculating  the  network  responses 
M  (x)  in  the  presence  of  a  parameter  change  was  developed  and 
Included  in  the  search  for  mrn||  M  (xj  -  M  ||.  This  is  described 
in  detail  in  Section  2.1 . 

(b)  Independent  Measurements.  Increasing  the  dimension  ni  of  M  (and 
also  of  MfxJ) by  adding  independent  measurements  should  improve 
the  accuracy  and  resolution  of  the  fault  analysis.  In  most 
practical  cases,  the  number  of  accessible  terminals  is  fixed,  so 
the  only  way  to  increase  the  number  of  measurements  is  to  use 

i 

more  than  one  frequency.  It  is  important  to  determine  if  the 
set  of  measurements  is  independent.  The  work  on  this  problem 
is  described  in  detail  in  Section  2.3. 

We  have  discussed  how  the  network  parameters  may  be  determined 
by  minimizing  the  norm  of  the  difference  between  the  measured 
and  computed  network  responses  | |  M  (x)  -  M  | | ,  where  M  and  M 
are  both  real  vectors  of  dimension  m.  In  the  discussion  which 
follows,  the  network  responses  are  considered  to  be  y-parameters 
at  the  accessible  terminals,  all  referred  to  a  common  (accessible) 
node.  Figure  1  shows  a  network  with  p  +  1  accessible  terminals 


from  which  we  define  the  pxp  y-parametcr  matrix 


y12 

y22 

V 


where,  as  usual , 


Vk  =  0,  k  f  r 


In  (3)  Is  and  Vr  are  the  phasor  response  and  source,  re- 

o 

spectively.  The  matrix  y_  therefore  represents  p  complex 


(2) 


(3) 


Figure  1  -  p-Port  Network  with  p+1  accessible  terminals. 

2  2 

parameters  or  2p  real  parameters.  These  2p  real  parameters 
form  the  vector  M  of  dimension  m,  i.e.,  m  =  dim  M  =  2p2,  an 
even  integer.  (The  exact  order  in  which  the  elements  of 
vector  M  are  formed  from  the  elements  of  matrix  is  not 

A 

important.  However,  the  measured  response  vector  M  is  formed 
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from  the  measured  y-parameter  matrix  £  by  the  same  rule  used  to 
form  one  computed  response  vector  M  from  the  computed  y-paraineter 
matrix  y) . 

We  close  this  section  by  displaying  a  formula  derived  by 
Johnson  [4],  [5],  which  will  be  used  in  subsequent  sections. 

Johnson  shows  that  the  pxp  matrix  of  measurable  y-parameters 
defined  above  is  given  by 

JL  =  ^  ^  [I  -  Aj  (A2Yb_A])-]  A^]  AT  (4) 

where 


is  the  node-to-branch  incidence  matrix  (as  defined  by 
and  Kuh  [9]  partitioned  into  the  accessible  (A^ )  and  the 
inaccessible  (A2)  nodes,  and  Yb  is  the  branch  admittance 
matrix  [4],  [5],  [10].  If  the  network  contains  b  branches, 

Yb  is  a  sparse  bxb  matrix  whose  elements  are  simple  functions 
of  the  branch  admittances  and  the  controlled  source  gains.  It 
therefore  plays  the  role  of  the  parameter  vector  x_  and  can  be 
expressed  in  terms  of 

2 . 1  Parameter  Changes 

In  order  to  find  min  ||  y_  (Yb)  -  y_  ||  it  is  important  to  compute 
lb 

y_  (Y b)  efficiently.  Such  an  algorithm  has  been  found  if  a  change  is 
made  in  only  one  parameter  at  a  time.  To  derive  this  result,  suppose 
Yb  is  changed  by  the  bxb  matrix  A  .  Then  from  (4)  the  y-parameter 
matrix  becomes 


1 
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I'  -  A}  (Yb+A)  tl  -  aJ  (Ag(YbfA)  aJ)"1  Ae(Yb+A)]  AT 


(!>) 


The  major  computational  difficulty  in  evaluating  the  right  side  of 
(5)  is  computing  the  inverse  of  the  new  inaccessible  nodal 
admittance  matrix  Y.>  =  A.,  (Yb+A)  l\l  .  Define  the  qxq  matrix 
Y.,  £  A.Y^  ,  where  q  is  the  number  of  inacessible  nodes.  Thus 
Y.J,  a  Yj,  +  A., A  and  we  wish  to  find  the  inverse  of  Y . 

During  this  investigation  we  have  discovered  that  if  only  one 
network  parameter  is  changed,  a  can  he  formed  as  the  outer  product 
of  two  vectors  u  £  col  (iij)  and  v  £  col  (v..),  i.e.  0  =  uyj  .  (In 
this  case  A^  has  unit  rank).  The  specific  way  in  which  u  and  v  are 
fonned  depend  on  whether  the  network  parameter  changed  is  a  passive 
admittance  or  a  controlled  source.  The  several  special  cases  are 
discussed  in  Sections  2.1.1  through  2.1.5,  which  the  reader  may  skip 
without  loss  of  continuity. 

The  product  A^a  may  therefore  he  expressed  as  the  outer 
product  ll_VT  ,  where  U  £  A,,u  and  V  =  A,,v  are  both  vectors  of  dimension 
q,  where  q  is  the  number  of  inaccessible  nodes.  We  may  now  write 


i'?_  *  y2  +  yy_  • 


(6) 


Householder  [11]  shows  that  the  inverse  of  Y^  can  be  formed  from 

Y^"1  ■  Y'1  [1  -  11(1  +  V^U)"1  V1  Y"1  j 

In  (7)  the  inverse  of  Y0  is  known  from  an  initial  computation  (in 
terms  of  the  nominal  parameter  values)  which  is  performed  only  once, 
and  1  +  V1  Y7,1  U  is  a  scaler.  Thus  the  right  side  of  (7)  may  be 
found  without  a  matrix  inversion. 


(7) 
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The  Fortran  code  for  computing  Y~^  us i ny  (7)  has  been  written 
and  may  be  found  in  Appendix  1  (SUBROUTINE  Y2PI2).  Sparse  matrix 
techniques  have  been  used  to  store  the  vectors  U  and  V. 

The  purpose  of  the  next  several  sections  (2.1.1  through  2.1.5 
is  to  show  that  a  change  in  any  network  parameter  changes  the  branch 
admittance  Yb  by  a  matrix  of  unit  rank,  i.e.,  where 

-  ”  y.  v.  *  anc*  li  s  col  (u.,),  y  t  C°1  (v..).  An  efficient  computer 
program  for  forming  the  vectors  u  and  v  using  sparse  matrix  techniques 
has  been  written,  and  the  listing  may  be  found  in  Appendix  1 
(SUBROUTINE  UVVEC). 

Johnson  and  Pennington  [10]  have  shown  that  in  a  network  con¬ 
taining  b  branches,  the  bxb  branch  admittance  matrix  can  be  formed 
from  the  element  matrices  as  follows: 

lb  ■  (1  +  l  \  -  l)']  (I  -  I  M  +  nn)  (8) 

where  1_  is  the  bxb  identity  matrix,  Y_  is  the  bxb  matrix  diag  (Y.) 
and  Y^  is  the  self-admittance  of  branch  i;  R  is  the  bxb  matrix 
(R.  .}  ,  where  R..  is  the  transistance  in  branch  j  controlled  by  (he 

'  J  •  J 

current  in  branch  i ;  «,  and  G  are  defined  similarly  and  contain 
branch  coupling  in  the  form  of  current  gains,  voltage  gains,  and 
trans conductances  respectively. 

2.1.1  Parameter  Change  is  a  Passive  Admit t one 

If  a  change  is  made  in  branch  admittance  k,  then  Y^ 
becomes  Y  +  d,  where  A  is  zero  except  for  clement  k.  A  little 
algebra  shows  that  Y^  has  been  changed  to  Y^  =  +  £,  where 


tThe  reader  may  omit  this  section  without  loss  of  continuity. 


A 


W  (AX)  -  (AU)  X  -  (AW)  (AX) 


(9) 


and  X  “  Y  -  Y  v  + 

W  =  1  -  Y  +  a 

AX  «  Ay  -  A  il 
AW  =  -A  R 

y  nit 

Equation  (9)  simplifies  since  (aU)  (aX)  =  o,  which  may  he 
shown  as  follows.  After  some  algebra,  using  (9), 

(AW)  (AX)  -  ( _AyRfl)  Ay)  (*.-  1), 

which  is  zero  if  Ay  R^  a  =  0.  Suppose  branch  admittance  k  is 

chanqed,  and  define  B  a  R  a  .  Therefore  13..  =  V  (R  ).  (a  )  .. 
3  -  hu  -y  i,i  a  v  in 7 1  m  v  y'mj 

But  all  elements  of  Ay  are  zero  except  for  element  (k,k). 


c  t)  ,  j  /  k 

Btk  ■  <*»>1k  <Ay>kk  •  v  f 


(1(1) 


Now  A ,R„  a„  =  A.,  B.  Consider  the  elements  of  column  k  of  this 
"jHii  -y  -y  — 

matrix: 


<4y 


Vrk 


(*y>rt 


ik 


(ID 


From  (10)  we  see  that  Ay  Rfl|  Ay  is  zero  except  for  elements  in 
column  k.  furthermore  the  summation  need  ho  carried  out  over 
element  i  =  r  only,  since  (Ay) r ■  e  0,  r  t  i.  Finally,  (a  )rr  E  0 
unless  r  °  k,  so  Ay  R)(  Ay  contains  only  one  non-zero  element, 

B«,  5y>kk  "  (Ay>kk(R1.,)kk(Ay>kk-  |!,,t  tMs  ,Ast 
expression  is  zero  because  branch  k  never  contains  a  trans- 

resistance  controlled  by  branch  k.  (This  is  modelled  instead 

as  a  simple  self  admittance).  QED. 
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_ ■■■■■  :<*r  H 


L 


After  deleting  the  last  term  in  (9),  the  change  in  the 


branch  admittance  matrix  Y  due  to  a  change  in  branch 
admittance  k  becomes 


A-Ay-AylL+aAy-aAyi. 


+  A  R  Y  +  A  R  C. 
-y  hii  -  ‘y  ~m  --m 


(12) 


*  A 

where  «  =  a  -  Y  R  is  the  Norton  equivalent  current 

*  A 

coupling  matrix  and  G  2  G  -Yu  is  the  Norton  equivalent 
transconductance  matrix  (see  Johnson  and  Pennington  [10]).  We 
will  now  show  that  all  terms  on  the  right  side  of  (12)  contribute 
elements  only  to  one  row  or  column  (except  the  term  p  _a  ,  q.v.), 
if  only  one  network  parameter  is  changed.  As  a  result,  A  may 
be  written  u  v\  as  explained  in  Section  2.1. 


1.  A  places  a  term  in  element  (k,k)  of  A. 

-y  “ 

2*  ^-y  ^ij  "  ^  ^Ay^im 

(Ay)ii  i'ij 

■  0  ,  i  /  k 

Therefore  Ay  u  places  elements  only  in  row  k  of  a, 


A  A 

^--^y^ij  "  jjj  aim  ^Ay 


=  ni.1  {Av)j.i 
-  0  ,  j  t  k 


Therefore  »  Ay  places  elements  in  column  k  of  a,  if  branch  k 
contains  current  controlled  sources. 

“  a1K  (^y  ^Kj 
"  aiK  (Ay*KK 
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If  such  a  term  were  present,  the  network  model  would  contain 
a  controlled  source  In  branch  i  ( a .  ^)  which  was  controlled  by 
the  branch  current  of  a  branch  containing  a  controlled  source 
(branch  k  contains  source  u^j).  Such  a  condition  is  not 
allowed  in  the  model  [10],  so  «  ]L  =  P  • 

The  final  two  terms  in  (12)  contain  the  product  a  R  . 

J  '“'111 

Element  (i ,  vi)of  this  matrix  is 
^ ij  ~  ^y^im 

■  <Vii  ■ 0  ■ f  * k 

liy  3*>X1  '  (Ay>kk  <WkJ  •  1  "  k 

We  now  consider  the  final  two  terms  in  (12). 

6-  <4y  5,,,  y,j  ■  l  <4,  »„ 

-inhj  Y,ij  «  0  ,  i  f  k 

h  ■  <vkk  (iwVj  vu  v  j 

Therefore  Y  contributes  terms  along  row  k  of  A  . 

6.  (iy  Em  ffchj  '  l  <iy  Vfn,  <Vmj  -°  •  '  '  k 

^■y  'm  ^  j  "  ~  ^-y  ^’nAij 

=  I  (\)kk  (,’nAm 
^Ay'kk  )•  ^mVj  Vj 


Therefore  A.  R  G  contributes  terms  alonq  row  k  of  a  , 

-y  -m  -m  — 

although  few,  if  any,  such  terms  arc  likely  to  exist  since 
the  presence  of  this  term  requires  that  the  network  contain 
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a  transresistance  controlled  by  the  current  in  a  branch  con¬ 
taining  a  voltaqe  controlled  source. 

2.1.2  Parameter  Change  is  ^  Trans  res  is tancc* 

If  a  change  is  made  in  transresistance  (R  then  R 

m  ij  -m 

becomes  R  +  A  ,  inhere  a  is  zero  except  for  element  (a  )... 
“in  1  ~f  r  ij 

Replacing  R)(  by  R^  +  A^  in  (8)  and  performing  some  algebraic 
reduction,  becomes  Y^  +  a  where 


A  =  -  Y  a  Y  -  Y  A  G 

—  - r  —  - r  -i:i 

Both  terms  in  (13)  contain  the  product  Y  A  . 


Expandi ng 


element  (r,  s)  of  Y  a^. 


tiir’rs  '  l  V™  <Ar>r 


Yrr  (Ar)rs  =  0  if  r  t  i ,  s  ^  j 


'V(i  <\><J 


Therefore  Y_  A(,  has  only  one  non-zero  element  (i ,  j). 
expand  element  (r,  s)  of  both  terms  of  (13). 


We  now 


>•  t*-*rX)rs  -l  O-Vrm  V 

•  <1  Vrs  ¥ss 


“  0  ,  r  f  i  ,  s  t  j 


tt-VDij-  <\>(j  Vjj 


t  The  reader  may  skip  this  section  without  loss  of  continuity. 


Therefore  Y  A  Y  contains  a  simile  non-zero  element  (i,  j). 
— r  — 


2-  <»  ir  ‘  I  <1  ir>™  >' V 


r  mi  in  ms 


0  ,  r  f  1 ,  m  f  j 


~  ^---r^ij  ^in^js 


’ll  <f’Jjs  •  V  S 


Therefore  Y  A  G  places  terms  in  row  i  of  A 

- r  — m  — 


2.1.3  Parameter  Change  is  a  Current  Gain 


If  a  change  is  made  in  current  gain  a . .  ,  then  a  is 

'  V ' 

changed  to  a  +  A  ,  where  A  is  zero  except,  for  element  (a  ).  • 
Replacing  n  by  «  +  a^  in  (8),  after  some  algebraic  reduction 
Y^  becomes  Yfa  +  A  ,  wliere 

A  =  -  A  Y  -  A  G 
—  —a  —  —a  -ill 

We  now  expand  element  (r,  s)  of  both  terms  of  (14). 

'  •  <A,  -V)rs  '  I  <\>rm  Vms 

’  <*„>rs  vss  ■  0  •  r  *  i>s  *  S 

^Ax  -hj  '  \jj 

Therefore  a  Y  contains  a  sinqle  non-zero  element  (i,  j)  . 

—a  — 

2.  (A*  Vrs  =  l  (Aa)rm  Vms 

m 

■  0  ,  r  M ,  m  f  j 

K  Vis'  Kh)  (Vjs  *  v  s 

Therefore  A  G  contributes  elements  in  row  i  of  A  . 

-xi  -m  ~ 


t  The  reader  may  skip  this  section  without  loss  of  continuity. 
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2.1.4  Parameter  Change  Is  a  Voltage  Gain* 

If  a  change  is  made  in  voltage  gain  (1 .  .  ,  then  £  is  changed 

•  J 

t0  H  +  »  where  A^  is  zero  except  for  element  (A^j.  Replacing 

£  by  £  +  a  in  (8),  after  some  algebraic  reduction  Y,  becomes 

H  0 

Y.  +  A  ,  where 

A=-Ya-«Ya  ( 

We  now  expand  element  (r,  s)  of  both  terms  of  (15). 

’•  <lVrs’JYnn«Vms 

=  Yrr  <Vrs  =  0  ’  r  *  1  •  s  *  J 

^ij  ~  Yii  K^ij 

Therefore  Y  a  contributes  only  one  term  to  A  . 

u  — 


2-  felVn-JV.'  *V« 


=  o  ,  in  M  ,  s  f  j 


^-^ri  ‘  «ri  (IVij 


=  ”r1  vii  •  v  r 

Therefore  £  Y.  A^  contributes  terms  down  columnj  of  4  . 

2.1.5  Parameter  Change  is  a_T ranscond^cj^ancef 

If  a  chanqe  is  made  in  transconductance  (G  then  G  is 

nrij  -m 

changed  to  G  +  /V  ,  where  A  is  zero  except  for  element  (a  ).  . 

-to  "9  "9  -g  i  j . 


t  The  reader  may  skip  this  section  without  loss  of  continuity. 


Replacing  by  0^  +  a  In  (b),  after  some  algebraic  reduction 


1  -  nn  J  -m  -g 

Y^  becomes  Y^  +  A  ,  where 


A  s  A  +  a  A 

-  -g  — g 


l\'e  now  discuss  both  terms  in  (lb) 


1.  A^  contributes  only  one  term,  (a^)^  ,  to  a  . 

2-  (»  4g»rs  *2  (»<,)„, s  =  0  .  m  i-  i.  s  /  ,i 

rtri  (4<ihj  ’  V  1 


Therefore  «  A^  contribute  terms  down  column  j  of  a  . 
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2.2  Finding  Independent  Measurements 

In  Section  1  we  discussed  the  problem  of  solving  the  nonlinear 

A 

equation  M(x)  «*  M  ,  in  which  x  is  the  unknown  parameter  vector  of 
dimension  n,  M  is  the  vector  of  m  network  responses  measured  at  the 
accessible  terminals,  and  M  (x)  is  the  vector  of  in  network  responses 
calculated  from  the  network  model.  The  three  cases  (m<n,  m=n ,  m>n) 
must  be  considered  separately.  Ransom  and  Saeks  [2]  have  discussed 
the  underdeterinined  case  (m<n),  which  is  not  very  practical  because 
the  fault  cannot  be  resolved  satisfactorily  in  this  case.  The  second 
and  third  cases  (m!n)  can  be  solved  by  finding  the  parameter  vector  x 
which  minimizes  the  norm  | |M(x)  -  M| |  .  (This  approach  has  the 
advantage  of  taking  measurement  error  into  account  in  a  sensible  way 
in  the  overdetemined  case,  m>n).  In  Section  2,  we  discussed  the 
solution  to  this  optimization  problem  by  means  of  a  search  over  the 
parameter  set  x_.  An  alternative  approach,  which  provides  some  insight 
into  the  question  of  independent  measurements,  is  to  expand  M(x)  in  a 
Taylor  series  about  the  nominal  parameter  set  .  If  we  truncate 
the  Taylor  series  after  the  linear  terms,  we  obtain  the  approximation 

M(x)  «  MtXg)  +  F(_x^)(x  -  xj 
where  Fjx)  is  the  Jacobian  Matrix  defined  by 

{F(x)l  « 

'  '  U  9Xj 

,  i  -  1 ,  2 ,  . . . ,  m 
,  j  =  1,  2 .  n  ( 

From  (16)  the  linearized  form  of  the  equation  M(x)  =  M  is  determined 

A 

to  be  F(xq)(x  -  x^)  -  M  -  M(Xq)  ,  in  which,  as  before,  dim  (x  -  x^)  =  n 
and  dim  (M-M)  =  m.  Since  we  are  treating  the  overdetermined  case  m>n, 


f  We  do  not  necessarily  recommend  using  (20)  to  determine  the  parameter  set  x 
of  the  faulted  network,  although  this  may  be  the  most  efficient  approach. 
See  Section  3  for  a  discussion  of  computational  complexity. 


-18- 


Since  min,  the  mxn  Jacobian  matrix  F  has  at  least  as  many  rows 
as  columns,  and  the  maximum  rank  of  F  is  n.  The  rank  of  F  is  less  than 
n  if  any  column  of  F  is  a  linear  combination  of  the  remaining  columns. 
The  simplest  case  is  that  a  column  is  a  constant  multiple  of  another 
column.  Suppose,  for  example,  (col  i)  =  a  (col  j),  i.e. 


ax^  a 


(k  =  1  ,  2,  . . . ,  m) 


One  obvious  way  this  can  occur  is  if  parameters  and  x.  are  two  like 
elements  in  series  or  parallel.  Suppose  that  x^  and  x.  are  the  con¬ 
ductances  of  two  resistors  in  parallel,  so  that  y  =  x>x.  15  the 
combined  conductance.  All  transfer  functions  which  do  not  access 
separately  the  inputs  or  outputs  of  x.  and  x.  nay  be  written  M,  (y,z), 

1  J  K  *  — 

where  z_  is  the  vector  of  all  network  parameters  except  x^  and  x^.  Now 

the  partial  derivatives  of  M,  with  respect  to  x.  and  x.  are 

x  1  J 


3fL  3fL  a ,  3M, 

_ n  _  _ _  __y  _ _ k 

3x^  3y"  3X.  “  ?y 


fV .  5. .  A  .  3 

3X.j  3.y  ax^  '  3y 

Thus  (21)  applies  with  a=l,  as  asserted.  This  result,  that  no  dis¬ 
tinction  can  be  made  between  imbedded  like  element  kinds  in  scries 
or  parallel,  is  intuitively  reasonable. 

A  less  obvious  example  of  column  dependency  is  illustrated  by 
the  simplified  cascaded  amplifier  circuit  shown  in  figure  2.  Define 
the  transfer  functions  and  network  parameters 


(21) 


(22) 

(2d) 


Figure  2  -  Simplified  Cascaded  Amplifier 


The  transimpedance  z,^  is  not  used  in  this  example  because  its  value 
is  zero.  The  Jacobian  matrix  (a'h/OXjl  is 


1  0  0 

9mlqm2R2R3  9ml9m2RlR3  9ml9m2RlR2 
0  0  1 


9m2RlR2R3 


9mlRlR2R3 


If  and  g^  are  equal,  then  columns  4  and  5  are  identical.  This 
doesn't  reduce  the  rank  of  £_  in  this  case  (the  rank  of  F  is  3),  but 
it  could  in  practical  cases  where  £  has  more  rows  than  columns  (m>n) 
Continuing  with  the  problem  of  column  dependency,  suppose  is 
the  conductance  of  a  resistor  in  parallel  with  a  capacitor  of 
capacitance  x.,  and  x  =  x.+jux.  is  the  combined  admittance.  All 

J  I  J 

transfer  functions  which  do  not  access  separately  the  inputs  or 
outputs  of  parameters  x.  and  x.  may  be  written  M.,(x,z),  where  z  is 

1  J  N  —  — 

the  vector  of  all  network  parameters  except  x,.  and  x^..  The  partial 
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i 


derivatives 

of  Mk 

with 

respect  to  x.. 

and 

x.  are 

J 

9Hk  . 

ax 

- 

3xi 

ax 

•  aXi 

ax 

DMk  _ 

d\_ 

ax 

aMk 

Jw  a'x " 

DX 

J 

ax 

•  ax  : 

Thus  (21)  applies 

wi  th 

a=jw  .  Columns  i 

and  j  of 

F  are  dependent. 

so  rank  £  < 

n  and 

the  nxn  matrix  F^F_ 

cannot  be  in. 

rted,  since 

rank  £TF  < 

n ,  if  n£n . 

As  a  consequence, 

(20)  cannot  be  solved  for 

the  new  parameter  vector.  We  have  discovered  the  surprising  result 
that  the  identification  problem  cannot  be  perfectly  resolved  if  the 
number  of  measurements  is  at  least  as  great  as  the  number  of 
parameters  and  the  network  contains  imbedded  elements  in  parallel 
(or  in  series).  That  rank  (£)  <  n  in  this  case  will  be  illustrated 
now  by  an  example. 

In  this  example  it  is  necessary  that  ,m>n  ,  and  we  have  chosen 
the  simple  2-port  network  illustrated  in  Figure  3  in  which  ni-n=3. 


Figure  3  -  Two-port  Network. 


We  choose  the  throe  transfer  functions  and  parameters 


n,  ■  ■  sr*  j»c  * r- 

fV,  "  y^j  B  JwC  +  6 

f'j  *  r  +  ® 

,  all.  ^ 

The  Jacobian  matrix  |  ^  is 

ju  1 

jut  1 

ju>  1 


0 

0 


x 


1 


c  l 


x3  "  0 


as  assorted,  columns  ?.  and  3  of  F  are  dependent .  Since  rank 
F  “  2,  the  fault  analysis  problem  cannot  be  fully  resolved  for 
this  problem. 

In  some  cases  the  column  dependency  can  be  removed  by 
increasing  the  number  of  measurements  by  using  another  frequency. 

In  th i s  example  if  Mj ,  and  are  defined  as  before  at  frequency 
w^w,  and  Mb  and  are  the  same  transfer  functions  at  frequency 
co=w.,,  the  Jacobian  matrix  is 


F  » 


7? 


jw^L 

0 

0 

1 

.irn^i 

o 

o 


j<*>i 


,iu' , 


j  w  ■; 


j"'2  ^ 

j<-'2  1 


and  rank  (F)  «  3. 


Another  way  to  handle  the  problem  is  to  split  the  transfer 
functions  into  their  real  and  imaginary  parts.  Continuing  the 
example,  define  the  transfer  functions  and  parameters 

=  Re(y-u)  =  G  x1  ■*  L 

M2  *  ^ll*  "  wC  -  i  x2  c  C 

M3  =  Ro(>'21^  G  x3  G 

a 

M4  c  Im(y21)  =  uC 

Now  m  =  4,  n  e  3,  and  the  Jacobian  matrix  is 

0  1 

(i)  0 

0  1 

w  0 

and  again  rank  (F)  c  3.  Since  and  are  identical,  rows  1  and  3 
of  F  are  the  same.  This  doesn't  reduce  the  rank  of  F,  but  it  is  clear 
that  may  be  discarded  without  reducing  the  information  of  the 
measurements. 

From  the  discussion  above  we  conclude  that  the  measurements 
should  be  split  into  their  real  and  imaginary  parts,  and  more  than 
one  frequency  should  be  used,  if  necessary,  to  make  the  number  of 
measurements  at  least  as  great  as  the  number  of  network  parameters 
(m>n). 

In  practical  examples  a  considerable  amount  of  computation  is 
needed  to  find  the  functional  form  of  the  elements  of  the  Jacobian 
matrix.  Sen  and  Saeks  use  such  functional  forms  in  their  discussion 
of  the  problem  of  independent  measurements  [6]  -  [0],  but  admit  that 
no  theory  yet  exists  for  choosing  measurements  while  taking  numerical 
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considerations  into  account. 

In  view  of  the  discussion  above,  it  appears  that  the  most 
practical  approach  to  the  problem  of  identifying  a  sufficient 
number  of  independent  measurements  (ni^n)  is  to  examine  the 

T 

numerical  conditioning  of  the  nxn  matrix  f  f  ,  evaluated  from 
the  nominal  network  parameters.  An  algorithm  for  evaluating  cotul 
(fTr)+  is  given  by  Torsythe  and  Moler  [12], 

Establishing  a  useful  set  of  independent  measurements 
requires  searching  for  min  cond  (fV)  over  various  sets  of  trial 
measurements  and  frequencies,  this  time  consuming  search  seems 
most  suitable  to  a  "simulation  before  test"  process;  i.e.,  the 
best  set  of  measurements  is  established  from  the  circuit  model 
before  testing  the  faulty  circuit.  The  contract  period  was  not 
sufficiently  long  to  permit  trying  this  procedure,  but  a  working 
algorithm  is  in  preparation. 


2.2.1  Computing  the  Jacobian  Matrix 

The  Jacobian  matrix  is  defined  by  f  ~  [aM. (x)/ax .] . 

We  have  taken  the  network  responses  H(x)  to  be  the  real  and 
imaginary  parts  of  the  y-parameters  at  the  accessible 
terminals,  where,  from  (A), 


*  •  MbA-i  -  W!:  feWi 


4  ■  U-iVi)  o:- 


(•' 


t  If  A  is  a  square  matrix,  cond  (A)  is  defined  to  he  the  ratio  of  tho  largest 
to  smallest  eigenvalue  in  absolute  value.  Therefore  cond  (A)  *  1,  and  large 
values  for  cond  (A)  indicate  nearly  singular  matrices  for  which  the  numerical 
evaluation  of  A-'  may  ho  meaningless. 
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and  so  we  want  to  compute  3y/ax^.  Making  use  of  the  identity 


5B"1  r  _  B-1  3B  (?-l 

3X  3X 

and  after  much  algebraic  simplification,  we  find  that 


-Vi  <WT>-'  U 


a!  (a^aV1  a. 


v }  a! 


(27) 


In  (27)  the  second  term  in  square  brackets, 

[1  -  A^  (A,,Y.  A..Y.  ]  ,  is  available  from  the  nominal 

—  "“'ll  — r-“~D — £  — C.—D 

analysis. t  The  first  term  in  square  brackets, 

Q  ”  Y(,A->  (A^Y^aJ)’  Ag]  ,  is  very  similar  to  the  second  term 
in  square  brackets.  In  fact,  if  the  network  model  contained 
no  controlled  sources,  Y^  would  be  symmetrical,  and  these  two 
terms  would  be  the  transpose  of  one  another.  The  right  side 
of  (27)  could  be  computed  very  efficiently  in  this  case, 
especially  since  oY^/dx^  is  very  sparse.  Unfortunately  we 
cannot  take  advantage  of  this  simplification,  because  the 
circuits  of  greatest  concern  contain  controlled  sources. 

An  algorithm  for  computing  the  right  side  of  (27)  as 
efficiently  as  possible,  taking  advantage  of  the  sparseness 
of  sY^/ax^,  has  been  written.  The  reader  is  referred  to  the 
listing  of  subroutine  UAf.COL  in  Appendix  1. 


^  The  term  appears  in  equation  (A)  and  is  computed  as  the  array  !)(*,•)  in 
subroutine  YMAT .  The  listing  may  be  found  in  Appendix  1. 
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3 .  Computational  Comp! cxi  ty 


Tho  algorithm  described  In  Section  2.1  makes  use  of  a  formula  of 
Householder  [11],  which  eliminates  the  need  to  invert  the  inaccessible 
nodal  admittance  matrix  Y'  to  calculate  the  accessible  transfer  functions 
^  each  time  a  parameter  is  changed.'*'  In  this  section  we  compare  the  amount 
of  computation  required  to  complete  the  search  for  min  ||M(x^)-H||  using 
this  efficient  algorithm  to  the  amount  of  computation  which  would  bo  re¬ 
quired  if  'T  were  directly  inverted  at  each  step  of  tho  search. 

The  computational  complexity  will  be  described  in  terms  of  the  number 
of  complex  multiplications  performed  by  the  algorithm.  This  depends  on 
the  network  configuration  and  specifically  on  the  number  of  accessible 
nodes,  inaccessible  nodes,  branches,  controlled  sources,  etc.  The 
computational  complexity  of  the  efficient  algorithm  is  even  more  difficult 
to  estimate  since  it  depends  in  addition,  on  such  measures  as  the  fraction 
of  branches  terminating  on  two  inaccessible  nodes,  the  fraction  of  branches 
terminating  on  only  one  inaccessible  node,  etc.  For  these  reasons  we  have 
made  the  following  assumptions:  the  total  number  of  network  parameters  n 
is  1.1  times  the  number  of  branches,  the  number  of  branches  is  2.2  times 
the  number  of  inaccessible  nodes,  AbX  of  the  branches  terminate  on  2 
inaccessible  nodes,  S5S  of  the  branches  terminate  on  exactly  1  inaccessible 
node.  Using  these  assumptions,  the  computational  complexity  can  be  calcu¬ 
lated  in  terms  of  tho  number  of  inaccessible  nodes,  q. 

The  basic  steps  required  to  compute  the  measurable  y-parameter  matrix 
£  are  the  same  for  both  the  direct  and  the  efficient  algorithms:  (1)  the 
new  branch  admittance  matrix  is  formed;  (2)  then  the  new  nodal  admittance 
matrix  Y^  is  formed  and  inverted;  (3)  finally  is  computed. 

*  See  equations  (4)  -  (7). 
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Steps  (1)  and  (3)  are  identical  in  both  the  efficient  and  the 
direct  methods.  Making  use  of  the  assumptions  discussed  above,  step 

9  . 

(1)  requires  about  O.Obq  complex  multiplications,  and  step  (3)  requires 

2  ? 
about  13.31q  complex  multiplications  for  a  total  of  about  1 3 . 3Gq  complex 

multiplications  for  steps  (1)  and  (3)  combined. 

Using  the  efficient  algorithm  to  invert  discussed  in  Section  2.1, 

step  (2)  requires  about  qL  +  3.4Gq  +  4.21  complex  multiplications.  Using 

the  direct  inversion  of  ,  step  (2)  requires  about  5q'  +  0.61q  complex 

multiplications. 

The  total  number  of  complex  multiplications  necessary  to  find  the 

measurable  y-parameters  y  after  a  single  parameter  is  changed  is  there- 
2 

fore  14.36q  +  3.4Gq  +  4.21  using  the  efficient  algorithm  and 
3  2 

5q  +  13.36q'  +  O.lGq  using  the  direct  algorithm.  The  efficient  algorithm 
is  clearly  very  superior. 

A  Fibonacci  search  using  7  evaluations  of  y  is  performed  for  each 
parameter,  and  the  number  of  parameters  is  assumed  to  be  1.1  times  the 
number  of  branches  or  about  2.42q.  Thus  the  number  of  complex  multipli¬ 
cations  required  to  search  over  all  parameters  is  16.94q  times  the  values 
cited  in  the  previous  paragraph.  These  results  are  shown  in  Figure  4  as 
a  function  of  the  number  of  inaccessible  nodes  q. 

We  close  this  section  with  a  discussion  of  the  computational 
complexity  of  forming  (20)  and  solving  for  the  unknown  parameters  x. 

Three  steps  are  required: 

(1)  form  the  Jacobian  matrix  F,  which  lias  a  column  for  each 
parameter  and  a  row  for  every  measurement; 

(2)  form  the  product  F.^F; 

(3)  solve  for  x  . 
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Using  the  assumptions  discussed  previously,  step  (1)  requires  about 
3  2 

9.29q  +  26. Mq  real  multiplications  per  parameter  for  a  total  of 
4  3 

22.48q  +  63.26q  real  multiplications.  Step  (2)  requires  about 

2.93q  m  real  multiplications,  where  m  is  the  number  of  measurements. 

Assuming  that  m  =  2n  :  4.84q,  then  step  (2)  requires  about  14.1 8q" 

real  multiplications,  finally  step  (3)  requires  about  4.72q'  +  5.86q  -  0.81q 

real  multiplications.  Therefore  the  total  number  of  complex  multiplications 

(1  complex  multiplication  =  4  real  multiplications)  to  form  (20)  and  solve 

4  3  2 

for  £  is  about  b.62q  +  20.54q  +  1.47q'.  This  is  also  plotted  in  figure  4. 

4.0  Conclusions 

Vic  have  studied  the  problems  of  high  resolution  fault  diagnosis  in  linear, 
analog  electrical  circuits  in  which  a  limited  number  of  terminals  are  accessible 
for  measurement.  He  have  chosen  to  describe  the  faults  in  terms  of  short  circuits 
between  pairs  of  nodes  and  the  values  of  parameters  in  the  network  model  (resis¬ 
tances,  capacitances,  inductances,  and  controlled  source  gains).  This  approach 
is  an  outgrowth  of  our  previous  study  [4],  [bj  in  which  the  fault,  description  was 
given  only  in  terms  of  short  circuits. 

Specifically  we  have  developed  and  tested  an  efficient  algorithm  for  computing 
network  responses  (transfer  functions  at  accessible  terminals)  after  a  single 
parameter  is  changed.  The  parameter  change  does  not  need  to  be  small,  and  in 
fact  the  algorithm  calculates  exactly  the  new  network  responses,  no  matter  how- 
large  the  parameter  change  is.  This  algorithm  was  developed  for  use  in  a  search 
over  the  network  parameters  for  the  network  responses  which  best  match  the 
measured  responses.  The  computational  complexity  (in  terms  of  the  number  of 
complex  multiplications)  of  this  algorithm,  which  takes  advantage  of  the  special 
relationship  between  the  network  before  and  after  a  parameter  is  changed,  was 
found  to  be  very  much  better  than  the  computational  complexity  of  a  direct 
algorithm.  (See  Figure  4). 
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In  addition  we  studied  the  problem  of  choosing  a  set  of  measurements  for 
fault  analysis.  We  have  found  that  the  suitability  of  a  set  of  measurements 
for  resolving  a  fault  can  be  expressed  in  terms  of  the  Jacobian  matrix  F  of  the 
network  responses  (the  partial  derivatives  of  the  network  responses  with  respect 
to  the  network  parameters).  The  fault  can  be  resolved  to  the  network  parameter 
level  only  if  the  number  of  measurements  m  is  at  least  as  .great  as  the  number  of 
network  parameters  n  and,  in  addition,  the  rank  of  r_  is  n.  This  latter  condition 
can  fail  in  the  various  ways  illustrated  in  Section  2.2.  The  most  promising  way 
to  find  a  set  of  independent  measurements  in  practical  cases  appears  to  be  to 
examine  the  numerical  conditioning  of  the  nxn  matrix  over  various  combinations 
of  measurements  (m  *  n).  We  have  written  an  efficient  algorithm  for  forming  the 
Jacobian  matrix.  Finding  the  "best-condi ticned" nxn  matrix  1-^1  appears  to  be  a  time 
consuming  process  which  seems  best  suited  to  a  "simulation-before-test"  process. 

While  investigating  the  measurement  problem  we  discovered  that  it  may  be 
possible  to  accurately  approximate  the  parameter  set  by  solving  a  set  of  linear 
equations  in  which  the  coefficient  matrix  is  fTF  .  The  computational  complexity 
of  forming  and  solving  these  equations  increases  as  (q  is  the  number  of  in¬ 
accessible  nodes),  whereas  the  computational  complexity  of  an  efficient  search 
3  ' 

increases  as  q  .  However,  Fiqure  4  shows  that  one  cycle  of  the  linear  approxima¬ 
tion  is  less  complex  than  one  cycle  of  an  efficient  search  (using  a  seven  point 
Fibonacci  search  in  each  coordinate  direction)  for  networks  having  fewer  than  40 
inaccessible  nodes.  Until  testing  is  completed,  however,  it  is  impossible  to 
judge  how  many  cycles  of  the  linear  approximation  are  needed  for  a  good  solution, 
and  problems  of  convergence  are  also  possible. 

5 . 0  Future  Work 

We  have  developed  several  efficient  algorithms  for  computing  network 
responses  and  partial  derivatives,  which  can  be  put  together  in  various  ways 
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to  solve  for  the  parameter  values  and  short  circuits,  if  any  exist.  Work  is 
currently  in  progress  to  evaluate  the  various  ways  these  algorithms  can  be 
used  in  forming  an  automatic  fault  diagnosis  tool  for  analog  electrical  circuits 
which  has  the  best  combination  of  efficiency  and  reliability.  This  process  will 
require  a  considerable  amount  of  testing,  including  the  use  of  real  measurements 
on  real  circuits. 
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"Efficient  Fault  Analysi 
in 

Linear  Analog  Circuits" 
by 

A.  T.  Johnson,  Jr. 
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Efficient  Fault  Analysis  in  Linear 
Analog  Circuits 

AU  lil  P  T.  .JOHNSON,  JR,  mimio  r.  ii  h 


Ahttmct- -I-'nuIi  In  hmhIoj’  nWw»»4i\  Is  u  ftwni  of  ik'ImoU* 

pM»nH't<*r  liVnlifKafion.  11k  pnjljlfiu  of  finding  {xmjuwcUM'v  (txnt) 

mcAXiirrmow*  *t  the  a«vc*\U>U'  terminals  can  Ik  cx|m'w.xi  hs  |)h'  v*a’u(nni 
of  a  xy stem  cf  lH»tllnci»r  rqtiatiwv  Niit'h  n  vwu  m  of  equation*  is  usually 
m4hxI  by  a  imiUMIttHHisKmal  souivN  t  xeix  s(v’;j  s\%  tiw  sv.uvb  w^nlrcs 
xohln#  for  (bo  ncfnwV  ix.'vp<v\vl'\  ht  low  of  (bo  |»a»an»r(»'»x,  o^iijviiiiij; 
(Ik  solution  to  (ho  OKflsorcsl  rrxponvs-s,  anil  fhow  (Ik  panPtiKtcx* 

in  such  a  wrx  as  to  prvxbuc  a  rcs^'nvc  closer  Ur  (Ik  tvc<ks\mM  rex)*on«c. 
The  computation  (Iok  required  (»>r  ibis  *nal)«U  |v  «»ficn  exvYwIvo  in 
practical  applications,  and  |n*se\  a  imxjoi  loyM>l»mcnl  to  fault  analyst*  In 
analog  circuits.  11k  represent  at  but  of  vlvnt  and  circuit  x  pens  add! 
Ilona  I  problems,  depending  upon  how  Ibe  iwtoori  rquatbrns  an'  for 
mutated. 

We  have  used  nodal  analysis,  tvjjrcsc«tltvj|  vlnni  circuits  b>  a  change  In 
the  nctmni  graph,  and  open  circuit*  by  n  rcio  (ranch  admltiarvY.  A 
formula  of  Householder  and  span*'  matrix  techniques  are  used  to 
efficiently  continue  fix*  response  oC  an  electrical  chvolt  oltb  either  a  sboit 
cfncuU  or  other  Uirc  pant  motet  change.  The  nuniJvr  of  complex  multi 
plica  lions  and  divisions  to  find  fl*e  netmuli  itsintn'c  at  ihe  accessible 
terminals  when  the  network  has  a  single  shon  circuit  Invoking  at  lem.t  one 
Inaccessible  terminal  is  reduced  by  shoot  a  facto*  of  IS  compnvd  with 
direct  methods. 

1.  INTRODUCTION 

T  is  surprising  perhaps  that  although  the  problem  of 
fault  isolation  in  analog  circuits  has  been  studied 
longer  than  the  comparable  problem  in  digital  circuits,' 
lire  digital  problem  is  much  better  understood,  and  auto¬ 
matic  test  equipment  (Al  l  )  for  isolating  faults  in  digital 
equipment  to  the  chip  and  even  to  the  gate  level  is 
available,  but  A  l  l-  available  for  isolating  faults  in  analog 
circuits  is  much  less  sophisticated  (1.1),  |bl|  t  he  explana¬ 
tion  of  this  phenominon  lies  partly  with  the  fact  that 
analog  signals  are  inherently  more  complex  than  digital 
signals.  They  occur  continuously  in  time,  rather  than  at 
discrete  times,  and  (heir  values  (in  principle!  have  infinite 
resolution,  instead  of  being  tinneakd  to  a  finite  number 
of  bits. 

In  the  linear  analog  circuits  discussed  in  tins  report,  the 
problem  of  signal  variety  can  I  e  partly  overcome  by 
restricting  obseivation  to  the  sinusoidal  steads  state,  since 
the  system  behavior  to  any  input  can  be  determined  fiotn 
a  knowledge  of  the  sinusoidal  steady  state  icsponsc  at  all 

M«mi*cupl  rfstnrd  StplrmN-t  I'D!*.  reviv'd  M.itvl'  S,  I 'IN  Hun 
work  s'  an  svipp»'ilrd  t'N  the  Air  I  otvr  Office  of  Scientific  RescArch 
I  hr.' Ugh  ihr  ISM  AS1  I  Nmmnr.  I  Acuity  KescAtvh  I'ropiAin  (V\  I’M  til 
under  Conlracl  1  ■»  ls.'lV7r>  ('.(Vin;,  I  lie  Ohio  Slate  I  Imeinn ,  Col 
umhm,  tilt 

The  author  n  with  llie  Onler  of  Tiumrcnna,  Wutener  i'ollrae, 
I'hfsier,  PA  l»UI » 

’An  estenvve  hil'hogrephi  of  l<olh  diRilal  and  analog  fault  isolation  is 
given  t»y  Sacks  and  I  itwily  (l/| 


frequencies.  Clearly  >i II  fieqneneies  cannot  be  considered, 
and  the  problem  of  how  to  choose  a  set  of  frequencies  has 
only  begun  to  receive  attention  |‘J).  |I0|  In  aslditixm.  the 
response  at  any  Irequcncy  is  relatcsi  iv«  the  paraiuctci 
values  and  the  netwoik  graph  by  tumlineai  equations. 
which  the  o»!\  method  of  solution  possible  is  a  time 
consuming  search 

In  this  paper  we  show  how  the  search  for  single  slum 
circuits  can  be  accomplished  efficiently  enough  to  make 
the  computation  practical.  We  also  dcmonstiatc  the 
effectiveness  of  this  method  in  identifying  short  circuit 
faults  in  analog  ci.cmis  whose  resistance,  capacitance, 
inductance,  and  controlle>l  source  gains  differ  si&mfi 
candy  from  their  nominal  values. 

II.  1  OUMUt  AttNCI  till  )  I’xRXMIMRSM  |||i 

Aix  t  ssii'ii  Nodi  s 

We  will  eonsidci  Imeai  lumpcdpaiameki  netwoik \ 
having/)  t  I  accessible  tciminals.  numbered  0  />,  as  shown 
in  Tig  I  The  network  of  Tie.  I  contains  n  4  I  nodes, 
including  ihc  reference  nods',  ansi  !'  branches  and  is 
assumed  is'  contain  ik>  internal  iiislepciisleitt  sources  but 
may  contain  depensleut  sources  We  will  formulate  the 
nodal  equations  using  ihe  standard  branch  shown  in  log 
a 

The  passive  clement  in  branch  a.  denoted  by  its  admit 
tancc  l'M.  must  be  nonzero,  I'tii  any  of  ihe  st'iiivs'  max  be 
zero  s'r  all  may  be  present,  if  desired  In  the  equations  that 
follow,  r  is  the  branch  voltage  veclsir,  j  o  tb,c  branch 
current  vector,  c  is  the  mvle  voltage  vector,  is  die 
independent  voltage  source  vector,  and  /  in  the  mdeprn 
stem  current  ssiuvee  \ectsir,  using  the  m'tatu'n  o(  I^cn.'ci 
and  Kuh  (1)  The  dimension  sif  each  s'f  these  svels'is  in 
hxl,  exsept  e  which  is  n  \  I  Considering  the  kill  stan 
dard  branch  shown  in  l-ig.  1,  it  is  clear  that  the  branch 
r  i  constraints  arc  given  by 

y »  f  plc  t  (o  YRm)j  >  »',  (1) 

where  >'^-sliag  (>,)  in  the  branch  pasoxe  ashmttanec 
matrix,  and  <i„.  /•,  o,  and  N„  ate  the  coupling  mainces 
due  to  Ihe  controlled  sources  The  dimension  of  each  of 
these  matrices  is  l'\b.  (Ihe  elements  of  matt  tees  and 
n  aie  curivnt  ssuirces,  ami  the  elements  of  malth-es  /i  and 
ate  voltage  somx'es  controlled.  rcs|vviisely,  by  soil¬ 
ages  and  cunents.)  l  ollowing  the  c»'mention  of  Desoer 
and  Kuh,  define  the  nx  I  innle  to  bianch  tnctilence 
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Fig.  1,  Linear  lumped-pararaetct  network  with  p  r- 1  accessible  termi¬ 
nals. 


Fig.  2.  Standard  branch  k. 


matrix  A, 

+  I.  if  branch  j  leaves  node  i 

{/<},,=  —I,  if  branchy  enters  node  i 

0,  if  branch  j  is  not  incident  on  node  i 

Solving  (I)  for  j  and  using  the  Ktrchhoff  current  law 
Aj*=0  and  the  Ktrchhoff  voltage  law  ,1  re  - t\  it  can  be 
shown  that  the  network  node  equations  are 

AYbA  re  =  >  (2) 

where  Yb  is  the  />  X  ft  branch  admittance  matrix  given  by 


branch  which  contains  a  current  controlled  source  docs 
not  control  any  source,  then  the  inverse  of  1+  Y'N^~a  is 
I  —  YK„  Hr.  Since  most  networks  can  lx*  modeled  this 
way,  forming  the  node  equations  dix's  not  present  any 
difficulty. 

The  accessible  nodes  have  been  numbered  0 -p,  and  the 
inaccessible  nodes  have  been  numbered  /.*  ♦  1  through  n. 
Partition  the  rows  of  the  nodc-to  branch  incidence  matrix 
A  into  accessible  and  inaccessible  nodes 


where  is  the  pX  n  incidence  matrix  of  accessible  nodes, 
and  A,  is  (he  qXn  incidence  matrix  of  inaccessible  nodes, 
and  q~n  -p  is  the  number  of  inaccessible  nodes.  Using 
(5)  in  (2)  and  partitioning  i  and  e  as  follows: 


where  i,  and  c,  arc  />  X  I  vectors,  and  i,  and  e,  are  q  X  I 
vectors,  we  obtain 

■W/'i  .'1,^/  lre, 

i  ,  .  (<>) 

j  {'I 

Now  eliminate  e,  from  (ft)  and  make  use  of  the  fact  that 
ij**0.  since  no  independent  sources  arc  connected  to 
inaccessible  nodes,  to  obtain 

'a2Y,\. 4  fa  (7) 

so  that 

'.42n]^,r  (S) 

is  the pXp  matrix  of  measurable /-parameters  defined  at 
the  ports  having  a  common  “ground"  node  (node  0  in  rig. 
1).  The  subscript  ;  indicates  that  the  /-parameters  arc 
evaluated  at  frequency  to,. 

III.  Fault  Analysis 

l.et  j.\  be  the  matrix  of  measured  /-parameters,  where 
the  measurement  is  made  at  frequency  to,.  A  fault  tcpic- 
sented  by  a  shift  of  the  network  element  values  from  their 
nominal  values  could  be  identified  by  solving  the  nonlin¬ 
ear  equation 


n-(l+  YRm~a)  \Y-Y,HGJ  (3) 

and  »  is  the  n  X  l  vector  of  independent  source  currents 

i-A(U-YRm-a)-'(Yv,~j,).  (4) 

Note  that  Vt\  is  a  Norton  transformation  of  the  indepen¬ 
dent  voltage  sources,  because  of  the  choice  of  standard 
branch  shown  in  big.  2.  the  network  can  always  be 
modeled  with  current  sources  only.. 

In  both  (3)  and  (4)  it  is  necessary  to  invert  the  tnutiix 
I  +yk„*  in  order  to  formulate  the  node  equations. 
Johnson  and  Pennington  |2)  have  shown  that  if  the 
network  can  he  modeled  so  that  the  branch  current  jk  of  a 


Yh)  (9) 

for  the  unknown  parameter  values.  If  a  single  frequency  is 
used,  as  in  (4),  the  number  of  network  parameters  usually 
exceeds  the  number  of  independent  v-parameters,  so  that 
(9)  has  no  unique  solution.  One  approach,  suggested  by 
Ransom  and  Kicks  |3|,  is  to  find  the  solution  tv’  (4)  that 
minimizes  a  norm  II  where  is  evaluated  from 

the  nominal  network  element  values.  Ransom  and  Sacks 
give  an  approximate  solution  to  this  problem  using  linear 
methods,  the  difficulty  with  this  approach  is  that  the 
solution  assumes,  roughly  speaking,  that  the  most  likely 
state  of  the  element  values  in  the  network  is  the  one  which 
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causes  the  smallest  drift  from  the  nominal  values  con¬ 
sistent  with  the  measurements.  This  assumption  excludes 
the  possibility  of  catastrophic  faults  (open  and  short 
circuits). 

Another  approach  is  to  augment  (9)  by  measurements 
at  a  sufficient  number  of  frequencies  that  a  unique  solu¬ 
tion  exists.  The  relationship  between  fault  resolution  and 
the  number  of  frequencies  and  measurements  is  discussed 
by  Sen  and  Sacks  [9J,  1 10).  Although  their  algorithm 
requires  choosing  a  number  of  frequencies,  it  is  not  clear 
how  this  is  done.  Thus  we  might  define: 

.v”(.v'i.»v  ■  •.»’*) 7  0°) 

and  solve  the  equation 

>’“>’(>  a).  (II) 

Nonlinear  equations,  such  as  (11),  are  often  solved  by 
defining  a  norm  c<™  ||  y-y ||,  and  then  making  a  multidi¬ 
mensional  search  for  min  e.  The  amount  of  computation 
required  by  such  a  search  could  be  prohibitive.  Chen  and 
Saeks  |4]  show  that  if  only  one  parameter  is  changed  at  a 
time,  the  measurable  system  function  can  be  evaluated 
without  the  need  to  invert  a  matrix,  which  is  apparently 
necessary  in  (8).  Their  method,  is  based  on  a  formula 
given  by  Householder  [5]. 

In  this  paper  we  show  how  the  Householder  formula 
can  be  used  to  efficiently  analyze  networks  containing 
short  circuits  by  considering  the  change  in  the  network 
graph.  The  measurable  effects  of  open  circuits  and  finite 
parameter  changes  can  be  analyzed  by  considering 
changes  in  Yh,  and  although  this  is  essentially  equivalent 
to  the  method  described  by  Chen  and  Saeks,  it  is  given 
here  in  the  context  of  nodal  analysis  so  that  all  paramelet 
changes  can  be  analyzed  with  one  formulation  of  the 
network  equations. 

A.  Short  Circuits 

If  two  nodes  of  a  network  become  shorted,  the  new 
node-to-branch  incidence  matrix  A'  is  related  to  the  origi¬ 
nal  incidence  matrix  in  a  very  simple  way.  If  nodes  i  and  j 
are  shorted,  then  A'  is  formed  by  replacing  row  i  by  the 
sums  of  rows  i  and  j,  and  then  deleting  row  j.  If  node  i  is 
the  ground  node,  then  A'  is  formed  by  deleting  row  j. 
These  row  transformations  can  be  obtained  by  prcmulti- 
plying  A  by  a  matrix  K ,  each  clement  of  which  is  either  0 
or  1.  Fig.  3  illustrates  the  row  transformation  matrix 
describing  two  different  short  circuits  (one  involving  the 
ground  node  and  the  other  not  involving  the  ground 
node)  in  a  network  with  5  nodes  and  7  branches.  The 
incidence  matrices  arc  partitioned  by  the  dashed  lines 
under  (he  assumption  that  nodes  3  and  4  in  the  original 
network  arc  inaccessible.  Thus  the  lop  half  of  the  A 
matrix  is  A ,  and  the  bottom  half  is  A}.  A  similar  partition 
is  shown  for  A’’m(A\  j  A])1,  but  in  the  transformed 
network  some  of  the  formerly  inaccessible  nodes  may 
have  been  made  accessible  (by  virtue  of  their  having  been 
shorted  to  an  Accessible  node). 
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pig.  3  Pffrd  of  shoit  circuits  on  nodc-lo  branch  iiKu’cme  moms.  A. 
(a)  Otiyinsl  network  ptaph.  (b)  Short  circuit  I'Ctwcrn  nodes  2  amt  .*. 
(c)  Short  circuit  between  nodes  2  and  ground  or  0. 


The  R  matrix  is  the  “aftcr-to-bcforc"  node  transforma¬ 
tion  matrix,  which  has  a  column  for  every  original  node  in 
the  network  and  a  row  for  every  node  after  the  fault. 
(Since  node  0  docs  not  enter  into  A,  it  docs  not  appeal  in 
R  cither.)  Thus  R  may  be  partitioned  by  columns  and 
rows  into  the  accessible  and  inaccessible  nodes  befotc 
(columns)  and  after  (rows)  the  fault,  as  shown  in  Fig.  3. 

Denote  the  partitioning  of  N  by 


R 


*11 

0 


C?) 


where  the  lower,  left  partition  in  (12)  is  0  since  the  faulted 
network  cannot  have  any  inaccessible  nodes  which  wcie 
accessible  in  the  unfaulicd  network.  In  other  words,  if  a 
short  circuit  occurs  between  an  accessible  and  an  inacccs- 
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sible  node,  the  combined  node  becomes  accessible.  Thus 

«■  RA  “ - -  - - 

*22^2 

and,  therefore, 

«*  ,  +  /? (y<2 

(13) 

B.  Forming  (he  y- Parameters  of  (he  Shorted  Network 

The  ^’-parameters  of  networks  containing  a  short  circuit 
can  be  computed  from  (S)  after  replacing/!,  by  A\  and  A2 
by  A\.  Define  Yi“AiYl>Af,  which  is  the  qXq  nodal 
admittance  matrix  of  inaccessible  nodes.  In  the  analysis 
that  follows,  we  assume  that  i\  is  known.  (In  fact  > j"  1 
must  be  computed  only  once  for  any  circuit  using  the 
nominal  clement  values,  so  the  cost  of  its  computation  is  a 
negligible  part  of  the  total  cost  of  analysis.) 

If  the  short  circuit  occurs  between  two  accessible  nodes, 
then  A'2  -A],  and  the  new  /-parameters  arc  found  from 
(8)  with  no  matrix  inversion.  In  fact,  if  a  short  circuit 
exists  between  two  accessible  nodes,  it  can  be  easily 
observed  by  making  resistance  measurements  between  all 
pairs  of  accessible  nodes,  so  the  method  described  here  is 
not  needed. 

If  the  short  circuit  involves  any  inaccessible  node,  then 
A2^A2  and  the  inverse  of  Y2-A2YhA'^  must  be  com¬ 
puted.  The  main  contribution  of  (his  paper  is  to  show  that 
the  inverse  of  Y2  may  be  found  in  terms  of  Y2  with  about 
one  fifteenth  ns  many  multiplications  as  are  required  to 
find  the  inverse  of  a  general  qXq  matrix.  This  saving, 
described  in  detail  in  the  next  section,  is  due  to  the  fact 
that  I  j  is  closely  related  to  >',,  whose  inverse  is  kown.  An 
additional  saving  of  computation  and  memory  is  realized 
by  the  use  of  sparse  matrix  techniques,  which  are  dis¬ 
cussed  in  Section  IV. 


C.  Finding  the  Inverse  of  YJ  Efficiently 
From  the  definitions  of  A2,  R22,  V3,  and  Vj  it  is  easy  to 
show  that 


y;~r22y2r[2- 


(14) 


Let  us  add  a  row  to  R22  as  follows 


S 


(15) 


where  uT  is  a  row  unit  vector  containing  all  zeros  except 
for  a  I  in  the  right-most  column  representing  a  shorted 
node.  (At  most  two  columns  of  R22  arc  involved  in  the 
short  circuit.  If  only  one  is  involved,  then  that  column  of 
u  r  is  unity.  If  two  columns  of  H22  arc  involved,  then  only 
the  right-most  one  is  unity.  This  convention  could  be 
reversed  by  changing  the  definition  of  B  )  We  now  decom¬ 
pose  S  into  V  and 


,V~t/+P 


(16) 


where  P  is  formed  ns  follows.  If  the  short  circuit  is 


between  an  accessible  and  an  inaccessible  node,  then 
/*«=■  G.  If  the  short  circuit  is  between  two  inaccessible 
nodes,  then  (P},_rJ_p  is  unity,  and  the  remaining  ele¬ 
ments  are  zero.  (In  this  case,  the  unit  element  in  P  is 
clement  { Af } , . .  where  node’s  i  and  j  arc  shorted,  and  i<j 
due  to  our  numbering  convention.)  With  these  definitions, 
it  is  easy  to  see  that  U  is  orthogonal,  since  its  columns  arc 
independent  unit  vectors  and  thus  form  an  orthonormal 
sc'  [7J.  Consequently  U  l/T,  and  if  the  short  circuit  is 
between  an  accessible  and  an  inaccessible  node,  then 
If  the  short  circuit  is  between  two  inaccessible 
nodes,  then  ,V  1  is  only  slightly  more  difficult  to  compute, 
as  wt  will  show.  Now  define  C  ■=  SY2SJ ,  and,  therefore, 


« 22' 

^'2 [  *22  j  ■»]- 

„  T 

U 

>  2  S  *22  >> 


UTY}R 22  :  uty2u 


and 


c  ~ 1 = (s r)  ~ 1  >  y  " 1  “ 


A’  i 

L 

a/T 

"n 

(17) 


Faddeeva  |8]  shows  that  A  =  >2"  1  [  1  —  (R22 F,u)A/].  and, 
therefore, 


Yi~  1  -  A[t~(7?j2)»Af]  '•  (18) 


The  inverse  on  the  right  side  of  (18)  can  be  computed 
efficiently  using  a  formula  given  by  Householder  |5],  |6). 
which  states  that  if  A  —  B  —  tnc 1 ,  where  v  and  n>  are 
column  vectors,  then 


A~'<=1 1 


B  W  rP>  1 
I  -  wrB  'v 


(19) 


where  the  division  on  the  right  side  of  (19)  is  by  a  scalar. 
Thus  the  inverse  of  Y2  is 


Yi~' 


(/?22  Y-Jl)M 
f  A/(A’,.J>) 


(~0) 


In  order  to  use  (20),  A  and  M  must  be  found  from  (17). 
We  have  mentioned  that  if  the  short  circuit  is  between  an 
accessible  node  and  an  inaccessible  node,  then  5  1  =  ,V 7 . 
In  fact,  if  the  short  circuit  is  between  an  accessible  node 
and  node  then  S-l.  and  K  and  A/  arc  simply  parti¬ 
tions  of  (17).  If  the  short  circuit  is  between  two  inaccessi¬ 
ble  nodes,  we  have  shown  that  S—U+P ,  where  17 
(tT,  and  the  rank  of  P  is  unity.  Thus  P  can  be  written  in 
the  form  t’ier,  and  (19)  can  be  used  to  find  S  '.  Proceed- 
ing  in  this  manner  it  can  lie  shown  that  S'  1  “  Ur  +  Q, 
where 


if  r s-n-p 
l  0,  otherwise 

And  so,  even  in  the  most  complex  case.  A"  and  M  can  be 
evaluated  without  a  matrix  inversion.  Some  time  can  be 
saved  in  computing  the  triple  matrix  product  in  (17)  by 
noting  that  /.  and  N  (each  column  vectors)  arc  not 
needed,  so  the  last  column  of  C  1  need  not  be  evaluated. 
Using  straightforward  matrix  multiplication  to  find  A'  and 


A 
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M  from  (17)  and  then  Y2  1  from  (20)  requires  q\q*-2qi 
+  3</J-</+l)  complex  multiplications,  where  q  is  the 
number  of  accessible  nodes  in  the  network. 

This  method  of  computing  ' 1  is  unnecessarily  in¬ 
efficient  since  many  of  the  matrices  in  (17)  and  (20)  are 
sparse,  and  others  (5  and  R22),  in  addition,  are  special 
since  their  non/.cro  elements  are  either  +1  or  —1.  In 
Secton  IV-A  we  show  how  to  reduce  the  number  of 
complex  multiplications  required  to  evaluate  Y2  1  to 
Pq*(q~  1),  where  ft  is  typically  about  0.3  and  is  always  a 
positive  number  less  than  unity. 

Once  Yj  1  is  known,  the  ^-parameters  at  the  accessible, 
nodes  are  found  from  (8).  Evaluating  the  matrix  products 
in  (8)  in  a  straightforward  manner  requires  2b*  +  (2p  +  q- 
\)b'  +  (q-l?b  multiplications,  where  b  is  the  number  of 
branches  in  the  network,  />  is  the  number  of  accessible 
nodes  not  counting  the  reference  node,  and  q  is  the 
number  of  inaccessible  nodes.  Using  the  sparse  matrix 
techniques  described  in  Section  IV-A,  the  number  of 
multiplications  needed  to  evaluate  y\  from  Y2  ~ 1  can  be 
reduced  to  2(1  +  y)b2,  where  y  is  typically  about  0.1.  For 
example,  a  network  with  20  branches,  3  accessible  nodes 
(plus  the  accessible  reference  node),  and  8  inaccessible 
nodes  requires  14180  multiplications  to  evaluate  y\  directly 
from  (8).  By  using  sparse  matrix  techniques,  this  is  re¬ 
duced  to  880  multiplications,  a  saving  of  a  factor  of  16  in 
this  case.  The  straightforward  evaluation  in  general  re¬ 
quites  nearly  h  times  as  many  multiplications  as  the  more 
efficient  approach  using  sparse  matrix  techniques,  so  the 
saving  increases  with  network  complexity. 

D.  Open  Circuits  and  Other  large  Parameter  Changes 

For  completeness,  we  include  a  discussion  of  open 
circuits  and  other  large  parameter  changes,  although  this 
has  been  discussed  elsewhere  [4],  The  main  purpose  in  this 
section  is  to  show  how  the  analysis  can  be  carried  out  in 
the  context  of  nodal  analysis. 

Parameter  changes  arc  accounted  for  in  (8)  by  a  change 
in  the  bxb  branch  admittance  matrix  Yh.  In  order  to 
evaluate  the  new y-parameters  from  (8),  the  inverse  of  the 
inaccessible  nodal  admittance  matrix  Y2  —  A }  YbA /  must 
be  found.  If  a  single  parameter  in  the  network  is  changed, 
then  the  new  branch  admittance  matrix  is  “  Yh  +  A, 
where  A  has  only  one  nonzero  element  representing  the 
change  in  one  parameter.  Suppose  a  change  occurs  in 
clement  ij.  Then  A  can  be  written  duv ',  where  d  is  a  scalar 
equal  to  the  parameter  change,  u  is  a  unit  vector  in 
direction  i,  and  v  is  a  unit  vector  in  direction  j.  Now 
define  U  “  A2u  and  V  &  A2v,  so  that  T2**  Yj  +  dllV1. 
Householder  shows  that  the  inverse  of  Y2  can  be  found 
from 

Yl-'-Y2-'[\-U(d-'+  VTYi'li)'  Vry,"'].  (21) 

The  inverse  of  l’2  is  known,  and  d~*+VrYjXV  is  a 
scalar,  so  (21)  can  be  evaluated  without  a  matrix  inver¬ 
sion. 


IV.  Sparse  Matrix  Techniques 

Sparse  matrices  arise  in  both  the  evaluation  of  the 
inverse  if  Y2  and  in  the  evaluation  of  the  accessible 
p-parametcr  matrix  yr  These  arc  discussed  separately  in 
Sections  IV-A  and  1V-B. 

A.  Evaluating  the  Inverse  of  Y2 

We  have  shown  that  the  inverse  of  Y2  may  be  found 
from  (20),  which  is  repeated  below  for  convenience: 

Yr '  -  A’l  >  +  ( >».”/(  1  -  >») ] 

where  K  and  M  are  defined  by 


The  dimensions  of  K  and  flf  are,  respectively,  (<7  —  1 )  x 
(q—  1)  and  1  X(q-  1). 

In  (2?)  Y{ 1  is  known,  and  S  is  defined  by 

(23) 

H 

R}]  is  the  matrix  which  transforms  A2,  the  (inaccessible 
node)-to-branch  incidence  matrix  before  the  fault,  into  A2. 
the  (inaccessible  nodc)-lo-brancli  incidence  matrix  after 
the  fault,  so  that 

d'2-R»A2.  (24) 

In  the  discussion  that  follows,  the  reader  may  find  it 
helpful  to  refer  to  the  illustration  of  Fig.  3.  Suppose  nodes 
i  and  j  are  shorted,  where  j  >  i.  Then  rows  1  through  i  -  1 
of  R22  contain  zeros  except  that  the  diagonal  element  is 
unity.  This  indicates  that  nodes  1  through  i--l  in  the 
faulted  network  arc  identical  to  nodes  1  through  »'— I  in 
the  unfaulted  network.  Row  i  of  R22  contains  all  zeros 
except  for  unity  elements  in  columns  i  and  j,  indicating 
that  ncale  i  in  the  faulted  network  is  formed  by  combining 
nodes  i  and  j  in  the  unfaulted  nctwoik.  Finally,  the 
remaining  rows  of  R22  contain  all  zeros,  except  that  the 
first  row  to  the  right  of  the  diagonal  is  unity.  This  indi¬ 
cates  that  the  remaining  nodes  of  the  faulted  network  are 
identical  to  the  remaining  nodes  of  the  unfaulted  network. 
The  faulted  network  contains  one  less  node  than  the 
unfaultcd  network,  so  tire  dimension  of  R22  is  (</  ~  1)  X  <7, 
where,  as  usual,  q  is  the  number  of  inaccessible  nodes. 

We  have  chosen  to  number  the  nodes  consecutively, 
beginning  with  the  accessible  nodes.  We  have  also 
avoided  analyzing  networks  with  a  short  circuit  between 
two  accessible  nodes,  since  such  a  fault  may  be  found 
from  driving  point  impedance  measurements  between  all 
pairs  of  accessible  nodes.  Thus  node  j  is  always  inaccessi¬ 
ble,  but  node  1  may  be  either  accessible  or  inaccessible.  In 
the  description  above  we  have  assumed  that  node  1  is 
inaccessible.  If  it  is  accessible,  then  is  simpler  than 
described  above,  since  all  rows  contain  a  single  unit 
element.  (The  unit  element  is  on  the  diagonal  in  rows  1 
through  j  -  I  and  in  the  first  column  to  the  right  of  the 
diagonal  in  all  other  rows.  This  indicates  that  all  inaccessi¬ 
ble  nodes  in  the  faulted  network  correspond  to  an  inacccs- 
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Fig.  4.  S  matrices  tor  a  network  with  8  accessible  nodes,  (a)  Short  Fig.  5.  (a)  Inverse  of  the  matrix  of  Tig.  -Ha),  (b)  Inverse  of  the  matrix 
circuit  between  nodes  4  and  b.  (b)  Short  circuit  between  nodes  6  and  of  Fig.  "Kb), 

any  accessible  node. 


sible  node  in  the  unfaultcd  network,  but  one  inaccessible 
node  in  the  unfaulted  network  is  now  accessible,  and  so  it 
has  no  row  in  R22.) 

In  (23)  u 1  is  a  row  vector  containing  all  zeros  except 
that  column  j  is  unity. 

To  illustrate  the  formation  of  the  S  matrix,  consider  a 
network  with  8  inaccessible  nodes.  (For  convenience  in 
this  discussion  they  are  numbered  1-8,  instead  of  />-)  1 
through  p  +  8.)  Fig.  4(a)  shows  the  S  matrix  for  a  short 
circuit  between  nodes  4  and  6,  and  Fig.  4(b)  shows  the  S 
matrix  for  a  short  circuit  between  node  6  and  any  accessi¬ 
ble  node.  The  formation  of  the  inverse  of  the  S  matrix  is 
given  by  5  " 1  «=  f / 7 .-+  Q ,  where 

fOl  =/  -*•  if  r™i~P'  s~n-  p 
^  ”  \  0,  otherwise 

as  explained  in  Section  111-C,  where  S*U+P,  and  all 
elements  of  P  are  zero  except  if  the  shorted  nodes  i  and  j 
are  both  inaccessible,  in  which  ease  [/’),,  is  unity.  The 
unity  element  of  P  is  circled  in  Fig.  4(a).  In  this  illustra¬ 
tion  we  have  taken  />“  0  for  convenience,  where  p  is  the 
number  of  accessible  nodes.  Also  note  that  n  -  p**q.  The 
inverse  of  the  S  matrix  can  be  formed  by  inspection,  and 
in  particular.  Fig.  5  shows  the  inverses  of  the  matrices  of 
Fig.  4. 

We  now  turn  to  the  problem  of  evaluating  K  and  Af. 
The  right  side  of  (22)  may  be  written 

( 1/+  Qt)  Y{  ’( UT+  Q) -  VY{  7/7  +  VY{ 'Q 

+  QrY{'UT+QrY2'Q.  (25) 


Recall  that  the  right-most  column  (column  </)  of  (25)  is 
not  needed.  The  second  and  fourth  terms  on  the  right  side 
of  (25)  are  zero  except  for  the  right-most  column,  so  they 
need  not  be  evaluated. 

To  evaluate  l/Y2  1  i / 7 ,  define  the  row  vector  V1  *=(t\|, 
where  vk  “(unit  column  of  row  k  in  U),  k  “  1,2,-  •  •  ,q~  1. 
Thus 

[Al  -[l/yf'l/H  “f  Vj"1],. 

l  J  m*r  L  *  J  mn  L  2  J  «„  ^  s  j  2***  O _  1 


Finally  Af  is  the  1  X(q  1)  row  vector  whose  elements 
arc  the  first  q  -  1  columns  of  the  last  row  (row  q)  or  the 
sum  of  terms  one  and  three  of  (25).  Thus 

~[  ],,*-[  >7' ],v  n-1,2,  ',*-1. 


A'  and  Af  ate  evaluated  from  (26)  and  (27)  with  no 
multiplications,  and  the  evaluation  of  A,  the  larger  of  the 
two.  requires  no  additions. 

Now  that  A  and  Af  are  known,  Y2~  1  can  be  found  from 
(21).  In  evaluating  (21)  it  is  useful  to  first  find  the  (q-  I) 
X  1  column  vector  jt  “  Ri2Y2u.  Recall  that  «r  is  a  vector 
with  all  zero  elements  except  that  row  j  is  unity.  Thus  x  is 
theyth  column  of  R}2  Y}.  R]}  is  the  partition  of  A  formed 
from  the  first  q  -  1  rows  (23).  We  may  make  use  of  the 
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row  vector  VT  to  form  x  ns  follows: 

[  Kj],,  r  if  node  i  is  accessible 

[-].-[*«»].-  [’)]*■'  ,°f  )• 

if  node  /  is  inaccessible 

«*=  1,2,-  •  •  ,q—  I.  (28) 

Using  the  definition  of  .t  in  (28),  (21)  may  be  written 
Vr,“A  +  A[vAf/(l -A/v)]  (29) 

which  requires  evaluating  three  matrix  products,  two  of 
which  directly  involve  x.  From  (28)  it  is  clear  that  x  is 
formed  by  reordering  tire  elements  of  column  j  of  T2. 
Since  only  a  few  inaccessible  nodes  arc  normally  coupled 
to  node  j,  x  is  usually  sparse.  Define  ft  equal  to  the 
fraction  of  nonzero  elements  of  column  j  of  K2.  Hie 
number  of  multiplications  required  to  compute  Mx  is 
fiq<Zq.  ft  is  usually  significantly  less  than  one,  and  so  it  is 
worth  designing  the  algorithm  to  avoid  the  multiplications 
by  zero.  Similarly,  pq  rows  of  xM  are  zero,  and  it  is  worth 
forming  only  the  nonzero  rows.  Thus  (29)  can  be 
evaluated  with  pq'(q-  1)  multiplications  and  divisions. 
Note  that  the  multiplications  in  (29)  arc.  the  only  multi¬ 
plications  required  to  evaluate  T2~  \  since  A'  and  A/  have 
been  formed  in  (26)  and  (27)  without  multiplications. 

B.  Evaluating  the  y- Parameters 

In  this  section  we  discuss  the  evaluation  of  the  y-param- 
eter  matrix  yt  given  by  (8).  In  the  faulted  network,  the 
node-to-branch  incidence  matrices  arc  A\  and  A'2,  so  that 
(8)  becomes 

»  =  /f  i  Y„[  \-A?Yf'A'2Yh}A?.  (30) 

Two  forms  appear  in  (30)  which  must  be  discussed  sep¬ 
arately:  A  tYA  and  A  YA  1 . 

The  most  convenient  way  to  store  the  nodc-to-branch 
incidence  matrices,  which  we  will  refer  to  generally  by  A 
in  the  discussion  which  follows,  is  by  two  column  vectors 
A*  and  A  “,  defined  as  follows.  If  the  Ath  column  of  A 
has  a  + 1  clement,  then  [A  +)t  ^  (row  number  of  the  +  1 
clement  of  the  Ath  column  of  A),  otherwise  [/f+]jk  =0. 
Similarly,  if  the  Ath  column  of  A  has  a  -  1  element,  then 
M“J*  “  (row  number  of  the  - 1  clement  of  the  Ath 
column  of  A),  otherwise  \A  ')*  *=0.  The  formation  of  A  + 
and  A  ~  corresponding  to  a  given  nodc-to-branch  inci¬ 
dence  matrix  A  is  illustrated  in  Fig.  6.  Note  that  in 
practice  it  is  necessary  to  store  only  the  vectors  A  *  and 
A  ~  and  not  the  matrix  A. 

Now  consider  the  evaluation  of  X^A  TYA .  Element 

I*L.,  is  g‘vcn  by 

A.  I 


-1  -1  1  0  0  0  0  0  0  c 

000  -1  100000 
0000011000 
00000-1-1  1  00 
C  0  0  0  -1  0  0  0  -1  1 


Fig.  6.  Node-to-branch  incidence  matrix  A,  ar.d  the  corresponding 
compact  equivalent  vectors  A  4  and  A  ~ . 


If  any  of  the  subscripts  of  ( T)  in  the  last  expression  is 
zero,  the  term  is  zero. 

The  evaluation  of  X'^AYA1  is  somewhat  different 
because  the  summation  is  along  the  rows  of  A  instead  of 
down  the  columns.  In  particular,  element  [A"]m  „  is  given 
by 

k  I 

“2MU-ML./[yk,-  <32> 

k.l 

The  right  side  of  (32)  cannot  be  evaluated  by  choosing  m 
and  n  as  in  (31).  Instead  A  and  /  are  chosen  and,  for  each 
choice,  in*.,  is  added  to  elements  and 

and  subtracted  from  elements  |A"]|^'|, .(.<  ], 
and  [XT^-j  i^.| .  If  any  of  the  subscripts  of  [A”]  is  zero, 
that  term  is  zero. 

Products  of  the  forms  X'  =  AYA7  and  X=ArYA  both 
appear  in  (30),  and  we  have  shown  how  to  evaluate  them 
without  multiplication.  An  additional  product  of  the  form 
X'  ■  Yb  must  be  evaluated,  where  Yh  is  given  by  (3).  Yb  is 
typically  very  sparse,  so  in  the  software  we  have  chosen  to 
stoic  it  as  three  vectors:  one  containing  the  nonzero 
elements  of  Yb,  and  the  other  two  containing,  respectively, 
the  row  and  column  of  that  clement.  Flic  saving  in  storage 
can  be  appreciated  by  recognizing  that  T„  is  a  6  x  h 
matrix,  where  b  is  the  number  of  branches  in  the  network. 
The  elements  of  Yb  which  are  not  zero  are  the  b  diagonal 
elements  plus  one  additional  element  for  each  controlled 
source.  If  the  number  of  controlled  sources  is  bf  10,  then 


t 


ir 


IMU!  TRANSACTIONS  ON  CIRCVITS  AND  SVS1TMS,  VO!.  CAS- 2b,  NO.  7,  JUI  T  1979 


Fig.  7.  Efficient  formation  of  the  matrix  product  X'-V,. 


the  sparse  storage  of  requires  storing  1.1  b  complex 
numbers  and  2.2b  integers,  compared  with  the  b2  complex 
numbers  required  for  conventional  storage.  If  the  software 
is  written  to  accommodate  networks  with  100  branches, 
then  conventional  storage  requires  setting  aside  room  for 
10*  complex  numbers,  compared  to  about  120  complex 
numbers  and  240  integers  required  by  compact  storage. 

Considerable  computation  time  is  also  saved  by  storing 
sparsely.  The  matrix  product  X‘  Yb  in  (30)  would 
require  bs  complex  multiplications  if  performed  directly, 
since  the  dimensions  of  both  A"  and  Yb  are  bxb.  The 
algorithm  of  Fig.  7  forms  the  product  in  (1  +  y)b2  complex 
multiplications,  where  y  is  the  fraction  of  the  b  branches 
containing  a  controlled  source.  Typically  y  is  about  0.1.  In 
Fig.  7,  YB(-)  is  a  vector  containing  the  1YB  nonzero 
elements  of  Yb,  and  YBROW(-)  and  YBCOl.(-)  arc  in¬ 
teger  vectors  containing  the  row  and  column  numbers  of 
each  clement  of  YB().  Tire  product  .Y'-V’,,  is  called 
/)(•,•).  D{- ,  •)  is  not  sparse,  so  it  is  stored  conventionally. 

Examination  of  (30)  reveals  another  matrix  product  of 
the  form  YpX,  where  in  this  case  Z=>  1 4-  A[Yj~  'A'2Yb. 
This  matrix  product  is  formed  in  a  similar  manner  to 
AT,  in  (l  +  y)b2  complex  multiplications.  Thus  altogether 
2(1  +  y)b2  complex  multiplications  arc  required  to  for¬ 
mulate  yr  A  direct  evaluation  of  (30)  requires  2b>-*-(q-  I 
+  p)b3  +  ((</-  l)J  +  /iJ)ft  multiplications.  Taking  />=»  2,  />” 
20,  1 1  wc  find  that  y,  requires  22880  multiplications  by 

direct  methods,  hut  only  880  multiplications  by  the 
methods  described  in  this  section  with  y»0.1,  a  saving  of 
a  factor  of  26  in  this  example. 

Before  closing  this  section  it  should  be  noted  that  the 
y-parameters  of  the  unfaultcd  network  are  evaluated  be¬ 
fore  evaluating  the  y-parameters  of  any  faulted  networks. 


This  is  done  by  applying  (30)  with  A\  replaced  by  A,.  A\ 
replaced  by  A2,  and  Yj  ~ 1  replaced  by  Yf  \  The  inverse  of 
Vj  must  be  computed  by  conventional  methods,  but  the 
methods  of  this  section  may  he  used  for  the  remaining 
computation. 

C.  Summary  of  Efficient  Computation 

The  y- parameters  at  the  accessible  ports  of  the  faulted 
network  arc  found  from  (30),  where  in  the  context  of  this 
report  the  faulted  network  is  formed  from  the  unfaultcd 
network  by  shorting  a  single  pair  of  nodes. 

In  Section  IV-A  we  showed  how  to  form  the  inverse  of 
the  (</—  l)X(<y-l)  matrix  l’J  with  only  (iq\q-\)  com¬ 
plex  multiplications,  where  /?  is  the  fraction  of  nonzero 
elements  in  Y2.  (Recall  that  Y2‘=A2YbA[  is  the  inaccessi¬ 
ble  nodal  admittance  matrix.)  Inversion  of  a  (<?—  1)M</- 
1)  matrix  requires  about  5(q  —  1)’  complex  multiplications 
using  NVilTs  rank  annihilation  method2  [11],  which  is  an 
efficient  application  of  the  Householder  formula  (19)  [5]- 

(7) 

In  Section  IV-B  we  showed  that  (30)  can  be  evaluated 
in  2(l  +  y)fcJ  complex  multiplications.  This  result  applies, 
whether  or  not  the  efficient  evaluation  of  the  inverse  of  l’j 
of  Section  IV-A  is  used.  Thus  the  total  number  of  com¬ 
plex  multiplications  required  to  evaluate  the  y-para meters 
at  the  accessible  ports  of  the  network  with  one  pair  of 
shorted  terminals  is  2(1  +  y )b2  +  q2(q~  1).  If  the  inverse  of 
Y2  is  found  by  conventional  methods  but  the  sparse 
matrix  techniques  of  Section  IV-B  are  used  to  form  y,, 
then  2(1  +  y)b2  +  5(q  —  1)J  complex  multiplications  are  re¬ 
quired. 

They-parametcrs  resulting  from  all  possible  single  short 
circuits  involving  at  least  one  inaccessible  node  must  be 
evaluated.  If  the  network  contains  />  accessible  terminals 
and  q  inaccessible  terminals,  then  it  can  be  shown  that 
A'  =  q(p  +  q)—  q(q  —  l)/2  distinct  pairs  of  shorted  termi¬ 
nals  can  be  identified.  The  total  number  of  complex 
multiplications  is,  therefore,  A7[2(  I  +  y)h2  +  Pq\q  -  1)] 
using  the  efficient  method  in  inverting  Yj  and  A'[2(l  + 
y)b2  +  b(q  —  l)5]  using  a  conventional  inversion  method. 
To  show  the  benefit  of  the  efficient  matrix  inversion,  let 
us  assume  that  the  number  of  accessible  ports  />  is  2,  the 
number  of  branches  b  is  twice  the  number  of  nodes, 
p  +  ^+l,  the  fraction  y  of  branches  having  controlled 
sources  is  0.1,  and  the  fraction  /3  of  nonzero  elements  of 
>'j  is  0.3.  Using  these  assumptions,  the  total  number  of 
complex  multiplications  necessary  to  analyze  all  the  in¬ 
accessible  short  circuits  using  both  the  efficient  and  the 
conventional  methods  of  inverting  are  shown  in  Fig.  8. 
The  advantage  of  the  efficient  method  is  especially  notica- 
ble  for  the  larger  networks. 

V.  Exampi.r 

Wc  have  described  an  efficient  way  of  computing  the 
y-parameters  at  the  accessible  terminals  of  a  network  with 

*An  trror  appears  in  our  edition  of  Ralston  »nd  Wilf  (ll|  In  the  flow 
churl  on  p.  75,  »iei  1 1  should  he  l-»i,  no!  2-«i. 
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Fig.  8.  Number  of  complex  multiplications  required  lo  anaiy7e  all 
inaccessible  short  circuits  in  a  network  with  q  inaccessible  nodes. 


an  internal  short  circuit.  The  question  remains  how  useful 
this  is  in  locating  short-circuit  faults,  since  the  network 
parameters  (resistance,  capacitance,  inductance,  and  con¬ 
trolled  source  gains)  of  the  faulted  network  will  be  differ¬ 
ent  from  their  nominal  values.  One  possible  way  to  locate 
the  fault  is  to  simulate  all  single  short  circuits  (involving  at 
least  one  inaccessible  terminal)  using  the  nominal  parame¬ 
ters  and  the  algorithm  described  in  this  paper,  compare 
each  of  the  simulated  y-parameters  to  the  measured  y- 
paranicters,  and  choose  the  closest  set  to  locate  the  short 
circuit.  Next  a  search  over  the  parameter  values  could  be 
conducted,  as  described  by  Chen  and  Sacks  [4]  or  in 
Section  Ill-D. 

To  date  the  short  circuit  algorithm  has  been  tested 
using  the  network  of  Fig.  9.  Measured  data  for  the  nine 
complex  y-parameters  at  the  three  accessible  ports  (nodes 
1-0,  2-0,  and  3-0)  were  obtained  by  simulating  single  shbrt 
circuits  and  setting  the  network  parameters  between  10 
and  20  percent  above  or  below  nominal.  It  was  expected 
that  this  might  completely  mask  the  effect  of  some  short 
circuits.  The  cost  function  used  to  compare  the 
'‘measured''  and  simulated  y-parameters  was 

v  k/- ->M 

f“£  no¬ 
where  y,j  is  the  y-parameter  simulated  from  the  fault 
model  and  y0  is  the  "measured"  y-parameter. 

Table  1  shows  the  actual  fault,  the  fault  determined  by 
the  minimum  cost  function,  and  the  value  of  the  cost 
function  of  the  correct  fault.  Of  the  seven  examples  at 
1000  Hr,  the  algorithm  correctly  located  five  pairs  of 
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Fig.  9.  Feedback  amplifier. 


TABLE  I 

Summary  or  Tesy  Ki.svi  ts  t 


shorted  nodes.  The  two  errors  scented  to  be  caused  by  the 
fact  that  the  magnitude  of  the  capacitor  impedances  at 
1000  Hi  were  much  lower  than  the  impedance  levels  of 
the  network  resistances.  (The  magnitude  of  the  impedance 
of  a  10-fiF  capacitor  at  1000  Hz  is  15.9  SI.)  Ii  was  expected 
that  the  short  circuit  between  nodes  6  and  7  would  lie 
difficult  it  identify,  since  the  50-jiF  bypass  capacitor  is 
only  about  3  12  at  1000  Hz,  and  so  the  changes  in  the 
parameters  of  the  measured  circuit  could  be  expected  to 
mask  the  short  circuit.  In  fact,  one  of  the  reasons  for  this 
test  was  to  determine  how  much  the  network  parameters 
could  be  changed  before  this  change  prevented  reliable 
identification  of  the  short  circuit.  In  this  example,  varia¬ 
tions  of  the  parameters  of  10  20  percent  from  nominal 
still  permitted  correct  identification  of  5  out  of  7  short 
circuits,  or  71  percent  of  the  eases  tried.  If  the  parameters 
of  the  “measured"  network  wete  set  to  the  nominal  val¬ 
ues,  then  correct  identification  would  occur,  of  course,  in 
100  percent  of  the  eases,  since  then  no  error  would  exist 
between  the  "measured”  and  simulated  y-parameters. 


pooooo 
pOuo oo 
poooc>o 


ooo  oooo^<h!!1  Li 

-*' o  oo o  O o  o o  O v >i  ’  ^ 

toot-h 


o  yo  0(0  OO  OO  O  ooo 


yOoo  cc-4 


Oo 


I  484 

VI.  Conclusions 

We  have  shown  an  efficient  way  to  analyze  networks 
containing  a  single  short  circuit  or  other  larger  parameter 
change.  The  efficiency  results  from  the  fact  that  the 
|  network  change  caused  by  cither  the  short  circuit  or  other 
large  parameter  change  results  in  an  alteration  of  the 
inaccessible  nodal  admittance  matrix  by  a  matrix  of  unit 
rank.  The  response  of  the  faulted  network  can  therefore 
be  computed  in  terms  of  the  nominal  response  with  many 
fewer  multiplications  than  are  required  for  and  original 
analysis.  Sparse  matrix  techniques  were  also  used  to  re¬ 
duce  execution  time. 

The  efficiency  of  the  algorithm  makes  it  practical  to 
exhaustively  search  all  possible  single  short  circuits  involv¬ 
ing  an  inaccessible  node  in  order  to  find  a  solution  to  the 
problem  of  locating  a  short  circuit  in  an  analog  network 
from  measurements  at  the  accessible  terminals.  The  test 


leKS  TRANSACTIONS  ON  CIRCUITS  AND  SYS1T.MS,  VOL.  CAS-it),  NO.  7.  JVl  V  1979 


(5)  A.  S.  Householder,  "A  survey  of  some  closed  methods  (or  invert¬ 
ing  lliRlnccs,"  J  SIAM.  vol.  V  pp.  I XX  Id 9.  |9X7. 

|0|  D.  Nohle,  Applied  /incite  Algebra.  1  nglcwood  Cliffs,  NJ.  then- 
tice-tlnll.  I9t.9,  pp.  147  148. 

[?}  - ,  Applied  Ijnear  Algebra.  Englewood  Cliffs,  NJ:  1‘ienlue- 

Hali.  1909.  p.  31.1. 

(8)  V.  N.  Ljiddeeva.  Computational  Methods  of  lent. 11  Algebra.  New 
York:  Dover.  1939,  pp.  102  103. 

(9)  N.  Sen  and  R.  Saeks,  "A  measure  of  testability  and  its  application 
to  test  point  selection- Theorv,"  t'r.o  loth  Muhcesi  ,Vi nir  Omuls 
Sys  Texas  Tech.  Umv.,  I  ubbock.  Texas,  Aug.  1977.  pp.  370  -383 

(|0)  - ,  **A  Measure  of  testability  and  its  application  to  test  point 

Selection-Computation,"  /‘roc  ACTOTFSTCObt  77.  pp.  212  219 
(Hyanius,  MA).  Nov.  1977. 

f  1 1  ]  A.  Ralston  and  It.  Wilf,  Miithenuiiica /  Meihi\ls  fiu  lhyil.il  t'om 
paters.  Vol  1  New  Yotk:  Wiley.  1907,  pp.  73-  77. 

[I2|  R.  Saeks  and  S.  R  Liberty,  Rational  Fault  Analysis.  Nets  York: 
Marcel  liekker.  |977. 

|I3]  ’T aigc  All:  sxsicms."  in  Cml  Fng.,  vol.  17,  no.  2,  March/April 
1978,  op.  28  -50 

[14J  E.  A.  lorrero.  "ATE-  Not  so  easy,"  1 FFF  Spectrum,  vol.  14.  pp. 
29-34,  Apr.  1977. 


examples  described  in  Section  V  show  that  short  circuits 
can  be  reliably  located  by  the  algorithm,  even  though 
measurements  were  made  at  only  one  frequency  and 
differences  between  the  actual  and  the  nominal  parameter 
values  were  not  taken  into  account.  We  expect  that  the 
reliability  can  be  improved  by  making  measurements  at 
more  than  one  frequency  and  by  solving  for  the  parameter 
values. 
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